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Abstract 



In this paper we present a natural and comprehensive generalisation of the stan- 
£> \ dard factorial moments (J-q) analysis of a multiplicity distribution. The Generalised 

Factorial Moments are defined for all q in the complex plane and, as far as the neg- 

■ ative part of its spectrum is concerned, could be useful for the study of infrared 
CO | structure of the Strong Interactions Theory of high energy interactions (LEP multi- 
plicity distribution under the Zq). The QCD calculation of the Generalised Factorial 

■ Moments for negative q is performed in the double leading log accuracy and is com- 
I pared to OPAL experimental data. The role played by the infrared cut-off of the 

Q H ' model is discussed and illustrated with a Monte Carlo calculation. 

(D 
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1 Introduction 

In this paper we present a natural generalisation of the standard factorial moments to 
continuous or fractional orders of a multiplicity distribution, which could be of interest in 
the study of the infrared structure of the strong interaction theory. 

In the past year, three groups of authors [I], fj, ^ have shown that it was possible 
to define and compute a multi fractal dimension, T> q , for QCD. Technically, this has 
been possible by computing the positive (and integer) order Factorial Moments of the 
distribution of particles in a restricted open angle A and could be compared with, say, 
the charged particle distribution in the Zq decay at LEPQ. 

J" (A) = ^(n-lUn-q + 1))^ ^ ^(i-^/d) (1) 
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where d is the dimension of the phase space under consideration (d = 2 for the whole 
angular phase space, and d = 1 if one has integrated over, say the azimutal angle). 

In the constant coupling case V q is well defined and reads: 

V q = 7o^— (2) 

where 7q = 4Cach s /2ti, a s is the strong interaction coupling constant, Ca is the gluon 
color factor. 

Dimensions T> q are called Fractal because they come from the natural generalisation 
to discrete variables of the standard moments which are used in the multifractal analysis 
of a continuous variable 0]. However, in this last field, the index q range is the whole real 
axis, while in our case it is restricted to positive integers. 

The choice of the factorial moments as a specific tool for the study of the scaling 
behaviour of the high energy multiplicity distributions have been of importance. As 
a matter of fact, it has been noticed by A.Bialas and R.Peschanski^ that the use of 
this observable permits to extract the dynamical signal from the Poisson noise in the 
Intermittency analysis of the multiplicity signal in high energy reactions. At first sight, 
the factorial moments T q are only defined on integer and positive values of q and does not 
gives any insight on the negative part of the multifractal spectrum (if any) of the nuclear 
matter. 

This has been noticed some years before by R.Hwa|7]] who first proposed a multifractal 
analysis of the signal by means of the so-called G-moments. However those moments did 
not have the property of the factorial moments to disentangle the Poisson noise from 
the dynamical signal, and thus suffer from statistical uncertainties. Further works are in 
progress in this direction 0. 

^From another point of view, one can understand, through their definition, that the 
standard Factorial Moments of the distribution are sensitive to the occurrence in the 
distribution of rare events of very high values of n as compared to its mean value For 
example, this is why the NA22 event [0] has been so important in the discovery of the 
intermittent properties of the high energy data. 

In contrast, the negative part of the q multifractal spectrum focusses itself on the 
study of rare events of relatively low values of the studied variable, which corresponds in 
our case to low multiplicity events (as compared to the mean value of the variable). The 
moments presented in the following have this property and are a natural and non trivial 
generalisation of the standard ones. 

However, in order to be efficient, one has to work with a multiplicity distribution with 
relatively high mean value, n « n&. This is why we will not apply this analysis to the 
intermittency analysis of the data but to the global multiplicity distribution and to its 
scaling properties with respect to the energy. 

Let us recall that the mean particle multiplicity produced by a gluon of energy E disin- 
tegrating in a cone of opening angle Bq is given, in the Double Leading-log Approximation 
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(DLA), by|nj: 

n b oc (3) 
were /i is the infrared cut-off of the theory, and that the corresponding global standard 



Factorial Moments follow by the KNOJTTJ phenomenon : 



F q = (n(n-l)..(n-q + l)) = c q n q b (4) 

where the Cq £1X6 known constants. At first sight, it could be difficult to 

understand that one could find out some scaling properties from moments which scales 
with rib- However, as it will be clear in the following, the Generalised Factorial Moments 
(GFM) analysis will show up a non trivial behaviour with energy. 

On the other hand it is known that the standard QCD factorial moments calculated in 
the DLA aproximation do not fit correctly the experimental data[|12|| and that important 



corrections are needed in order to describe the experiment reasonably . Some progress has 
been recently made in this direction Jl3|. In this paper we will however restrict ourselves 



to the DLA approximation of the theory for the calculation of the GFM of negative 
arguments of QCD and we will study the effect of the infrared cut-off of the theory on 
this observable. We will work both analytically and numerically with a simple QCD 
Monte-Carlo model (based on the Fragmentation structure [|I], |14], |I5|1). We also restrict 
ourself to the fixed coupling constant case. In all this calculation, we did not introduce 
any ad-hoc parameter exept the perturbative coupling constant, 70 which we fixe around 
.5 at the Z . 

As a consequence of these restrictions, we did not try to make any precise comparison 
with the experimental data. We leave to further works a more precise study of the sub- 
leading QCD corrections, running coupling constant and/or non perturbative corrections 
which may be important in this specific field. 

The paper is organised as follows : In section 2 we present the Generalised Factorial 
Moments (GFM) and their behaviour for some useful examples such as Poisson distribu- 
tion, self-similar structure (KNO) or Negative Binomial distributions. Then, in section 
3, we discuss the Generalized Factorial Moments of QCD in the Double Leading-log Ap- 
proximation. Then we discuss the introduction of the infrared cut-off in the theory. We 
conclude in section 4. 



2 The Generalised Factorial Moments. 
2.1 Definition. 

The standard Factorial Moments of a multiplicity distribution P n are given by 

F q = (n(n-l)(n-2)(...)(n-g + l)) P 
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= £p n n(7i-l)(...)(n-g + l) (5) 
o 

which, using the properties of the T (Euler) function can be writen as : 

and under this form can be continued in the complex q plane . 

The importance of the Factorial Moments comes mainly from the fact that they can 
be derived from the generating function G(z); from the theoretical point of view it is in 
general much easier to handle than the multiplicity itself. One has : 



g G(z) 
q dz« 



G(z) = J2z n P n (7) 







which one has to generalise to continuous or fractional values of q. 

Let us recall here the properties of the principal value distribution x^^ 1 defined on 
[0, 00 [ with respect to the convolution of functions defined on the positive real axis (causal 
functions) [l(J : 

T(-q) Jo r{-q) qJ[ ] [ } 

where / is any well behaved function, continuous and indefinitely differentiable at x=0. 
When the exponent q is positive, the action of x~ q ~ x (principal value) on a test function 
(f is given by : 



x q 1 , ip(x)) = I X 



- c '- 1 ^(x)-{^(0)+x^(0) + ... + ^^}), n g = lnt(q) (9) 



where Int(g) is the integer part of the real part of q. 

When the real part of q is negative, the integral is defined and can be calculated. 
When q is positive or 0, the principal value ansatz @ must be used, and, on positives 
integers the result is just the q th derivative of the function. This comes from the fact 
that for q integer and positive the integral diverges together with the T function at the 
denominator of [8[ 

Using this definition of the derivative one can define : 



d q G(z)\ z=i = YT^f Q G ( l - t ) t ~ q ~ ldt ( 10 ) 



In order to be complete, one has to prove that this definition is consistent with that of 
Eq. (|7|). It is easy to verify that introducing the definition of the generating function G 
in Eq. (P~0| ) one recovers Eq. (Q) thanks to the property of the (B) Euler function of the 
second kind. 
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2.2 Examples. 

i) The Poisson case. In this case, P n is given by : 



n , 



= — r exp ( - n 6 ), = exp (n b (z - 1)) 

n! 

where n& is the mean value of n, and it is easy to obtain : 



'11} 




a) b) 
Figure 1: The Poisson GFM; a) : n b —l, b) : n b =10 



T 



Fjnl 



J*(-q,n b ) 



(12) 



where j{— q, n b ) is the incomplete 7 function and 7* the analytical incomplete 7 fun- 
ction[[D]]. In any case, the value of 7* on the integer and positive (or 0) values of q is 1 
which is natural since one recovers here the standard factorial moments of the Poisson 
distribution. But one has to notice that, besides those points, the shape and behaviour 
of this function depends drastically on the n b value; This function has in fact two types 
of behaviour. One for n b ~ 1 and one for n b » 1. This is illustrated in figure |I[ 

If n b = 1 say, figure (|l| a) exhibits a steeply oscillating behaviour for q > and goes 
rapidly to when q is negative. This behaviour could prevent us from using those moments 
for small values of n b where one cannot wait for a faithful behaviour of the moments and 
where the numerical formula (|7|) can be very unstable. This is not a surprise if one 
considers that those moments are devoted to the study of rare events of low n, ie n « n b 
For n b » 1, say n b = 10, figure ([I] b) shows that the 7* function is practically 1 in a 
large interval of q : — n b « q « 3n b . 

^From the point of view of the intermittency data, this indicates that the Poisson 
noise will be disentangled from the dynamical signal only for those q greater than — n b . 
As in those data the mean value of the number of particle in each bin tends rapidly to 0, 
one can understand that the fractional part of the moments does not gives any dynamical 
insight on the basic process. 
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Figure 2: The GFM of the Negative Binomial Distribution, a) q behaviour for n& = 10; 
upper curve k — 1, intermediate one : k = 5, lower one : k — 10; b) rz^ behaviour at fixed 
q = — 5, upper curve fc = 1, intermediate one : k = 5, lower one : = 4 



zz ^ 5e// similar distributions. The best way to buid an asymptotic self-similar distribu- 
tion is to construct P n as a compound Poisson distribution (a particular case of these 
distributions is the NBD distribution). At the level of the generating function, this gives: 

roo 

G(z) = #0-1), H(u) = / (fix) exp (ux)dx, (13) 

J o 

where p is the usual KNO function : 

P n ~ lim ip(n/n b )/n b , (14) 

rij, —>oo 

and one gets: 

iJ(-u) = h(urib) 



: — A_ / hi-u^u-^du (15) 
T(-g) Jo 

On this integral, one can notice that if n b is sufficently large, T q = F q /n\ tends to a 
constant and we recover the KNO result provides h(u)u~ q — > when m — > oo. Unless this 
condition is fulfilled, jF g will depends on rib- 

If, say, h(—u) — u a , when w — > oo, JF 9 will be KNO for g >> —a. When g < —a, the 
integral in [I5| will be dominated by the high u behaviour of the integrand and 

T q ~ V + * (16) 

Measuring the generalised moments of the distribution provides a rather nice tool for the 
study of the high u behaviour of the generating function. 

As an exemple of this situation, we have calculated the moments of the Negative 
Binomial distribution : 

fc M = JT^W (17) 
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which gives : 



2 F 1 (k, -q, l-g, -n b /k) 
- q) 



where 2 -Pi is the hypergeometric function of the second kind|p~7|. 

The general trend of the reduced generalised moments of the distribution is given in 
figure (H-a) while the power- like behaviour of the moments is shown in fig (Q-b). 



3 QCD Generalised Moments. 

3.1 QCD Double Leading-log Approximation. 

The QCD evolution equation for the generating function G(Q, z) of the multiplicity dis- 
tribution produced by a parton of energy E disintegrating in a cone of opening Oq has 



been calculated since a while in the DLA approximation[10] and reads : 



dG{Q,z) r 1 2 {G(xQ,z) - l)dx 

~ G{Q,z) / 7 (19) 



dlog(Q) 'Jo ,u x 

where Q = EQq is the hardness scale. Notice that an implicit infrared cut-off must be 
understood in this equation, Qx > Qq. This cut-off tends to in the weak coupling regime 
of the equation and will be of importance in section 3-2. 

Let us first fix some notations. As in the preceding section, we define H(Q,u) = 
G(Q,1 — u), and the self similar solution (KNO) of the equation [19], h(v), such as 
H(Q,u) = h(ribu). Further, for negative u, let us define h(y) = h(—exp (y)). 

With these notations the QCD solution obeys an integro-differential equation which 
does not depends explicitely on the coupling constant 70, and reads (using ref.4 with some 
slight change of notations) : 

d 2 log(h) 

dy 2 

This equation has been solved in an implicit way (for negative u = — exp (y) |1 

rX(y) du 

y-y+ 



(20) 



/io g (2) ^2(u-l+exp(-tt)) 
X(y) = Ml/My)), y + = -.251 (21) 

Starting from formula (|T5|), the QCD GFM for negative q can be given as an y integral: 

1 /-logOib) 

This expression is well adapted to the steepest descent technique which gives: 

-X(y)-qy~-(y-y*)/2a 2 -l, 
y* = -q -c Q , a 2 = 1 + exp ( - g 2 /2 - 1), c = .41 (23) 



1 /-log(n b ) 

^ = W^)l-oo ^{-X{y)-qy)dy (22) 
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where cq has been numerically computed for y > 2. This gives asymptotically (n^ — > oo): 

log(^f d ) = g 2 /2 - (1 - c )g - 1 + (g - .5) log(-g) + log(<r 2 )/2, 

g « -2 (24) 

where we use the same approximation (steepest descent) for the T function. One has 
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Figure 3: GFM for negative values of q. Continuous line : Asymptotic QCD; dotted- 
dashed line : finite corrections to asymptotic theory; squares : finite infrared corrections 
(Monte-Carlo); rounds : OPAL data (one hemisphere multiplicity) 

now to deal with finite rib corrections to this formula. First we verify that the Gaussian 
approximation is a lower bound for h, the same for the the GFM calculated with formula 



22| . Using formula ^ with the Gaussian approximation for h(y) as a lower bound for the 
finite size QCD GFM gives : 

(logK) + c + q)/V2^] \ 



erf [ 



(25) 



The theoretical prediction are shown in figure |3| together with the experimental data from 
OPAL[|l!|. We have taken the multiplicity data on one hemisphere in order to be as close 
as possible to the one parton multiplicity distribution. Notice that the prediction of the 
theory including finite rif, effects has to be understood as a lower bound for the moments, 
specially for small values of —q. 



3.2 Infrared effects. 

Let us now turn back to formula (jT9|). We have now to take into account the finite infrared 
cut-off, Q of the theory. 
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In order to be coherent with the probabilistic interpretation of QCD DLA, the decay 
probability, U(x) of a parton of hardness Q into two partons of hardness xQ and (1 — x)Q, 
must be normalised positive and definite . Keeping tracks only for the logarithmically 
divergent part of the Altarelli Parisi kernel, we have : 



U(z) 



6(,) + £ 



lim(M 



z-e 



+ 8(z)(l 



7o lo £ 



'Me)) 



where Q(x) is the Heaviside step function. 

In the fixed (and finite) coupling constant case, one can observe that the probability 
distribution is no more positive definite in the limit e — > 0. As a consequance, this limit 
is incompatible with the building of a Monte-Carlo calculation. This is a reflection of the 
fact that the perturbative theory is exact only in the limit where 70 — > 0, i.e. at infinite 
energy. 

To evade this difficulty, one has to use a finite cut-off Q . The simplest choice is 
given by \og(Q/Qo) ~ I/70 wich gives II(z) = 7q@(<2 — exp (1/7q))/z. Using this form in 
equation (|l9| ) gives at the level of the self-similar solution h : 



OPAL 




Figure 4: n b behaviour of JF_ 5 : Continuous curve JF_ 5 oc nf, grey squares 
Monte-Carlo, dotted-dashed one QCD lower bound for finite size effects. 
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U(y) 
Hv) 

which has an asymptotic solution : 



y 1 

h(z)dz 

y- 1/70 7o 



(26) 



h(y) oc exp 



To' 



y » 



To 



(27) 
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in the limit where y goes to infinity. 

This exponential behaviour is different from the Gaussian solution of the QCD case 
(Effi). For finite 70, the infrared cut-off have changed the asymptotic behaviour of the 
self-similar solution of the QCD evolution equation. Notice however that one recovers the 
previous solution when 70 goes to 0. 

As the GFM of negative order are sensitive to the high y values of the generating 
function, it is not astonishing that they show a different energy behaviour : 

T q oc n~ q ~^ = (28) 

This behaviour is of the NBD type (see Eq. (|HJ)) is well reproduced by numerical (Monte- 
Carlo) calculations, even for low n&. 

Let us now discuss a little the QCD results. In figure [| the Monte-Carlo results are 
shown to be higher than the lower bound we derived in the weak coupling limit. We verify 
that this does not depends on the detail of the Monte-Carlo. The infrared structure of the 
theory, namely the cut-off in the evolution equation, enhances the fluctuation pattern for 
negative values of q. Notice that the Monte-Carlo calculation gives the known results for 
positive values of q which shows that the infrared cut-off has no particular effect on the 
standard positive moments. The OPAL data have been presented in this figure to show 
that the data (one hemisphere data) are substantifically lower than DLA QCD predictions 
in all the q range. Figure f| shows the n& behaviour of the GFM of order -5. It is interesting 
no notice that while the two QCD results are different in strength, the apparent slope of 
the QCD bound is close to the Monte-Carlo one and is very close to the one predicted in 
formula (H)(-g - I/70 ~ 3 at LEP) . 



4 Conclusions. 

In this paper we have presented a comprehensive generalisation of the standard factorial 
moments which have been shown to be sensitive to the infrared structure of the theory. 
This analysis focusses on the low multiplicity events of high energy reactions such as Z 
decay at LEP. The QCD GFM have been calculated together with a bound on low energy 
corrections and have been shown to be substancially higher than the OPAL data. When 
one includs the natural infrared cut-off in the theory, the asymptotic picture of the GFM 
are modified in the negative part of its spectrum and enhances the fluctuation pattern. 
Our feeling is that, in contradiction to the positive q case where Next to Leading-Log cor- 
rections (energy conservation effects [DJ) are probably enough to explain the discrepancy 



of the asymptotic theory with experiments, the GFM of negative order will emphasis the 
non perturbative part of the theory and could provide a glance on the hadronisation part 
of the strong interaction theory. Further work is in progress in this direction. 
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